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Evolution  of  Zinc  Oxide  Nanostructures  from  Non-Equilibrium  Deposition  Conditions 

PI:  Xudong  Wang,  Materials  Science  and  Engineering,  University  of  Wisconsin-Madison 


Introduction 

This  project  focuses  on  to  quantitatively  understand  the  time-dependent  evolution  of  Zinc 
Oxide  (ZnO)  nanostructures  driven  by  non-equilibrium  deposition  conditions,  and  thereby 
establish  an  experimentally  grounded  foundation  for  designing  nanostructures  and/or 
functionalities  via  kinetic  control.  Morphology  is  one  essential  element  that  gives  rise  to 
extraordinary  physical,  chemical,  and  mechanical  properties  in  nanomaterials.  Precise 
morphology  control  of  nanomaterials  is  a  notorious  task,  which  heavily  relies  on  fundamental 
understanding  of  the  governing  atomistic  mechanisms  and  kinetics  at  the  nanoscale.  Through 
vapor-solid  deposition  under  dynamically  controlled  conditions,  we  observed  three  unique  ID 
and  2D  nanocrystal  growth  behavior,  which  revealed  a  series  of  unique  growth  kinetics  at  the 
nanoscale  under  conditions  far  away  from  equilibrium.  First,  the  wedding  cake  growth 
mechanism  in  the  formation  of  ID  and  2D  ZnO  nanostructures  was  observed  within  a  narrow 
growth  window.  An  interesting  ID  to  2D  morphology  transition  was  also  found  during  the 
wedding  cake  growth,  when  the  adatoms  overcome  the  Ehrlich-Schwoebel  (ES)  barrier  along 
the  edge  of  the  top  crystal  facet  triggered  by  lowering  the  supersaturation.  Furthermore,  an 
inverted  wedding  cake  growth  phenomenon  was  observed  in  2D  nanostructure  evolution.  The 
atomic  steps  were  nucleated  along  the  edge  and  propagated  toward  the  center  forming  a  unique 
basin-shaped  crown  on  the  tip  of  each  nanowire.  The  morphology  and  temperature  relation 
revealed  the  ES  barrier  of  ZnO  to  be  0.88  eV.  Based  on  the  layer-by-layer  growth  model,  a 
unique  high-energy  facet  evolution  phenomenon  at  the  NW  tip  was  observed  at  different 
deposition  supersaturation  within  a  narrow  vapor  deposition  window.  The  facets  evolution  and 
exposure  of  high-energy  facets  were  attributed  to  the  fluctuation  of  the  energy  barriers  for  the 
formation  of  different  crystal  facets  during  the  growth  of  NW  tip.  This  project  substantially 
enriches  our  understanding  on  the  fundamental  kinetics  of  nanostructured  crystal  growth  and 
provides  a  transformative  strategy  to  achieve  rational  design  and  control  of  nanoscale  geometry. 


1.  Observation  of  wedding  cake  growth  mechanism 

We  reported  an  observation  of  the  wedding  cake  growth  mechanism  in  the  formation  of  ID 
and  2D  ZnO  nanostructures.  The  wedding  cake  growth  model  yields  a  very  similar  terraced 
surface  as  that  generated  by  a  screw  dislocation  but  does  not  possess  a  dislocation  core  and  spiral 
pattern.  An  interesting  ID  to  2D  morphology  transition  was  also  found  during  the  wedding  cake 
growth,  when  the  adatoms  overcome  the  Ehrlich-Schwoebel  (ES)  barrier  along  the  edge  of  the 
top  crystal  facet,  triggered  by  lowering  the  supersaturation. 


1 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Figure  1.  Morphology  of  microplate-caped  micropillars.  A.  An  individual  microplate-caped 
micropillar  clearly  shows  a  thread-like  feature  on  both  plate  and  pillar  surfaces.  B.  Top  view  of  a 
microplate  exhibiting  a  hexagonal  outline  with  multiple  concentric  terrace  circles.  C.  Tilted  and 
enlarged  view  of  a  microplate  showing  the  intermediate  sub-terraces  extruding  along  the  six 
equivalent  corner  directions.  D.  A  flat  and  smooth  central  terrace  surface  without  any  observable 
screw  dislocation  steps.  E.  Side  view  of  a  microplate  showing  the  terraced  structure  descending 
from  the  center  of  the  plate.  Inset  shows  the  flat  top  terrace.  F.  An  image  of  a  very  thin 
microplate  revealing  orientation  of  the  microplate  and  micropillar  is  differentiated  by  30°.  G.  and 
H.  Examples  of  plate-caped  micropillars  with  small  (G)  and  large  (H)  plate  thicknesses.  I. 
Irregular  terrace  patterns  due  to  the  shadowing  effect  at  high  plate  density  region. 

The  unique  ZnO  microstructure  was  obtained  under  a  dynamic  deposition  condition,  where 
the  growth  chamber  temperature  was  rapidly  increased  from  room  temperature  to  1673  K  within 
90  min.  During  the  ramp,  both  the  Zn  vapor  concentration  and  the  ZnO  deposition  temperature 
rose  continuously.  The  ZnO  microstructure  composes  a  micrometer- sized  pillar  that  was  capped 
with  a  large  hexagonal  microplate  (Fig.  1  A).  The  pillars  were  typically  several  micrometers  long 
with  a  diameter  of  ~1  micrometer.  The  microplates  were  ~2-3  micrometers  in  diameter  with 
variable  thicknesses.  A  thread-like  feature  can  be  clearly  observed  both  along  the  side  of  the 
pillar  and  on  the  plate  surface.  The  top  microplate  exhibited  a  sharp  hexagonal  outline  with 
multiple  concentric  terraces  (Fig.  IB).  These  terraces,  however,  did  not  exhibit  a  sharp 
hexagonal  outline  as  the  microplate  showed,  whereas  extrusions  together  with  intermediate  sub¬ 
terraces  (indication  of  rapid  growth  front)  can  be  observed  along  the  six  equivalent  corner 
directions  (Fig.  1C).  Although  this  surface  feature  is  very  much  similar  to  a  screw  dislocation 
growth  pattern,  a  spiral  center,  which  is  the  characteristic  feature  of  a  screw  dislocation-driven 
growth,  could  not  be  observed.  All  the  plates  exhibited  a  flat  and  smooth  central  terrace  surface 
(Fig.  ID).  This  phenomenon  fits  well  to  the  wedding  cake  growth  model,  where  new  atomic 
layers  nucleate  on  the  current  top  layer  giving  rise  to  a  wedding  cake  like  structure.  A  side  view 
of  the  plate  clearly  displays  the  step-wise  terrace  configuration  (Fig.  IE),  indicating  the  terraces 
were  evolved  initially  at  the  center  and  then  spread  out  to  the  edge.  A  zoomed-in  image  also 
evidences  the  flat  top  terrace  without  a  spiral  center  (inset  of  Fig.  IE). 


2 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


From  very  thin  plates  that  electrons  can  penetrate  through,  the  orientation  relationship 
between  the  plate  and  pillar  can  be  observed.  As  shown  in  Fig.  IF,  the  darker  hexagon  at  the 
center  of  the  plate  outlines  the  pillar  underneath.  The  hexagonal  capping  plate  exhibited  a  30° 
rotation  in  its  orientation.  Based  on  the  crystallography  of  wurtzite  ZnO,  this  feature  suggests 
that  the  six  equivalent  side  facets  of  the  pillar  and  the  plate  are  either  one  of  the  two  low-energy 
surfaces  of  {01  10}  and  {11  2  0},  respectively.  The  thickness  of  the  plates  exhibited  a  large 
variation  ranging  from  50  nm  to  ~  1  pm  (Fig.  1G,  H).  Regardless  of  the  thickness,  all  the  plates 
had  a  very  uniform  width  along  their  entire  thickness  with  a  clear  convex  terraced  top  surface. 
While  most  plates  exhibited  a  six-fold  symmetric  terrace  pattern,  irregular  terrace  patterns  could 
be  observed  at  high  density  region  (Fig.  II).  Shadowing  effect  or  competition  of  precursor  might 
be  responsible  for  the  asymmetric  growth  behavior. 

The  crystal  structure  of  ZnO  plate-caped  micropillars  was  characterized  by  Cs-corrected 
scanning  transmission  electron  microscopy  (STEM).  Selective  area  electron  diffraction  (SAED) 
pattern  taken  from  the  cross-sectional  structure  revealed  that  the  entire  structure  was  a  single 
crystal  and  grew  along  the  <0001>  direction.  The  height  of  each  terrace  was  in  a  range  of  5~20 
nm.  The  terrace  edge  was  not  perpendicular  to  the  plate  surface  and  all  of  them  showed  a  slope 
of  -153°.  High  resolution  STEM  (HR-STEM)  images  taken  at  the  plate-pillar  joint  region 
showed  a  coherent  crystal  lattice  and  no  grain  boundary  or  other  dislocations  could  be  found. 
Particularly,  no  indication  of  the  existent  of  a  screw  dislocation  core  was  discovered  along  the 
center  position  where  the  top  terraces  nucleated.  The  terrace  edge  was  found  composed  of 
multiple  small  steps  without  any  dislocations,  consistent  with  SEM  observations.  The  HRTEM 
observations  further  evidence  that  the  wedding  cake  growth  model  governed  the  evolution  of 
ZnO  crystals. 

Supported  by  microscopy  observations,  the  wedding  cake  growth  model  is  used  to  analyze 
the  evolution  of  the  2D  microplates.  In  a  wedding  cake  growth,  formation  of  surface  terraces 
relies  on  nucleation  of  new  atomic  layers  on  the  top  surface  due  to  the  accumulation  of 
concentrated  adatoms  confined  by  the  energy  barrier  at  step  edges.  This  yields  the  representative 
morphology  of  concentric  circular  mounds  without  a  spiral  feature.  The  wedding  cake  growth 
behavior  was  mostly  observed  in  the  low  temperature  epitaxial  growth  of  metal  surfaces. 
However,  in  our  hybrid  structure,  evolution  of  suspended  2D  plates  requires  the  lateral  growth  to 
exceed  the  supporting  pillar’s  top  surface.  Therefore,  the  adatoms  need  to  overcome  the  ES 
barrier  which  prevents  adatoms  from  diffusing  down  at  the  terrace  edge  and  then  incorporated  at 
the  corner  along  the  vertical  facets.  Lattice  kinetic  Monte  Carlo  (LKMC)  method  was  used  to 
simulate  the  morphology  evolution  when  the  ES  barrier  becomes  tangible.  This  simulation 
verifies  that  the  ES  barrier  triggers  the  formation  of  the  suspended  plate  morphology  through 
decreasing  the  supersaturation. 

To  better  understand  the  wedding  cake  growth  mechanism  in  nanostructure  evolution, 
numeric  calculations  were  implemented  to  analyze  the  growth  rates  along  different  directions  in 
comparison  to  the  screw  dislocation-driven  growth  mechanism.  In  a  wedding  cake  growth,  the 
lateral  growth  rate  v/  is  a  measure  of  how  fast  adatoms  attach  to  the  ledge  front,  which  is  same  as 
that  in  a  screw  dislocation  growth.  Quantitatively,  vi  can  be  presented  following  the  Burton- 
Cabrera-Frank  (BCF)  theory: 
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V,  = 


2/1  v  exp  (-^Zl) 
kT 

l  +  -eXp(2AU~^“) 


(1) 


where  a  is  the  distance  between  adsorption  sites,  v  is  the  vibration  frequency,  (pm  is  the  work  to 
detach  an  atom  from  half-crystal  position  to  the  vapor,  k  is  the  Boltzmann  constant,  T  is 
temperature,  <r  is  the  supersaturation,  As  is  the  mean  free  path  of  an  adatom,  8  is  average  spacing 
of  the  kink  sites,  AU  is  the  kinetic  barrier  for  the  incorporation  of  the  building  units  into  the  half 
crystal,  and  Edes  is  the  desorption  energy  from  the  surface  to  the  vapor. 

As  for  the  vertical  growth  rate  va,  in  the  wedding  cake  model,  it  is  directly  related  to  the 
nucleation  rate  of  new  layers,  which  occurs  only  when  the  top  layer  reaches  a  critical  size. 
Therefore,  the  vertical  growth  rate  can  be  formulated  as: 


v„ =  ■ 


h_ 

R 


(2) 


where  Rc  is  the  critical  size  of  the  top  layer  and  h  is  the  height  of  the  newly  formed  terrace.  On 
the  other  hand,  the  vertical  growth  rate  of  the  screw  dislocation  model  follows  the  BCF  theory: 


=  avexp(- 


^1/2 


kT  cr 


)  —  tanh- 


(3) 


where  crc  is  characteristic  supersaturation.  Since  there  is  no  nucleation  involved  in  the  screw 
dislocation  growth,  wedding  cake  model  usually  exhibits  a  slower  vertical  growth  rate. 

One  direct  result  of  the  different  axial  and  lateral  growth  rates  is  the  number  of  terraces  on 
the  plate  surface.  Shown  in  Figure  2A,  within  the  deposition  area  (~2  mm  in  length)  where 
terraced  plates  were  received,  the  number  of  terraces  decreased  as  the  sampling  point  moved 
along  the  gas  flow  (downstream)  direction,  whereas  no  obvious  size  variation  could  be  identified. 

Then  we  converted  the  number  of  terraces  into  growth  rate  ratio.  By  rearranging  the 
expression  of  lateral  growth  rate  (Eq.  1),  the  relation  is  displayed  as  below: 


vi  = 


2  ov  exp(-^^-) 
kT 

l  +  -eXp(2AU~^) 


X. 


(4) 


For  any  given  plate  structure,  we  assume  a  and  T  are  constant,  and  thus  the  lateral  growth 
rate  is  linearly  related  to  2S.  Therefore,  when  Xs  decreases,  v/  decreases. 

To  determine  axial  to  lateral  growth  rate  ratio,  we  assume  the  time  spent  by  the  first  terrace 
(the  base  terrace)  to  grow  from  the  plate  center  to  the  plate  edge  (//)  equals  to  the  time  to  evolve 
the  rest  terraces  (ta).  The  base  terrace  expands  with  a  velocity  of  v/j  and  then  decreases  to  17,2 
when  approaching  the  edge  due  to  the  decrease  of  effective  adatom  collection  area  (?/  =  t/j  +  t/j). 
The  height  from  the  base  terrace  to  the  plate  top  is  nh,  where  n  is  the  number  of  the  terraces  and 
h  is  the  height  of  one  step.  Thus, 
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nh  R-  / l 


v. 


'  t 


1,2 


vl,l 


(5) 


where  R  is  the  plate  radius.  From  the  equation  of  the  lateral  growth  rate  (Eq.  4),  we  can  get  the 
expression  for  v/,  ? : 


vu=A\ 


(6) 


where  the  coefficient  A  is  a  constant  by  assuming  that  a  and  T  are  constant  during  the  growth  of 
a  specific  microplate: 


A  = 


2ovexp(-^) 
_  kT_ 

,  5  2AU-L 

1+  exp( - 

a  2kT 


'des  ' 


As  for 

vi,2  =A(As-x) 

where  x  is  the  distance  the  terrace  grows  with  the  velocity  vv, 2 •  (2*  ~  x)  represents  the  distance 
from  the  terrace  edge  to  the  plate  edge.  This  equation  expresses  the  instantaneous  velocity  at  the 
point  (A s  -  x )  because  v/,2  varies  with  the  effective  collection  area,  v/,2  can  be  expressed  in  terms 
of  vij  as: 

A-x 

Vl,2=—7 - Vl,l 

As  (8) 


Then  t/>2  can  be  obtained  by  integrating  from  x  =  0  to  x  =  (As  -  d): 


x-d  dx  cx-d  dx  A „  ,  2 


_  rA-d  ax  _  r 

L2  “Jo  ~  “  Jo 

Vl,2 


=  ^-1n 

A-x  V,,  d 


2. 


V,  , 

-  V  M 

VU 


Insert  Eq.  9  into  Eq.  5,  we  get: 


nh  R-As  A  A 

—  = - -  +  —  In  — 

d 


vu  vu 


Rearrange  Eq.  10  to  get  the  growth  rate  ratio: 

va  _  nh 

v , 


A, 


"  (R-As)  +  As  In  — 
d 


(9) 


(10) 


(11) 


The  values  of  n,  h  and  d  can  be  obtained  by  measurement,  therefore  the  number  of  terraces  n 
is  converted  to  the  growth  rate  ratio  vjvij- 

Next,  the  position  was  converted  to  the  supersaturation.  The  local  vapor  pressure  is  obtained 
by  subtracting  the  amount  consumed  by  deposition  from  the  initial  amount  yielded  by  the 


5 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


precursor  source.  According  to  the  previous  results,  the  deposition  of  ZnO  from  the  vapor  phase 
underwent  two  steps,  zinc  condensation  and  zinc  oxidation.  At  high  temperature  (>700  K),  the 
condensation  of  zinc  is  the  rate-limiting  step.  Therefore,  we  use  the  condensation  rate  of  zinc  as 
the  deposition  rate.  The  condensation  rate  is  determined  by  the  sticking  coefficient  (S)  and  the 
incident  molecular  flux  (F)  via  Rcon  =  FxS,  where  S  is  estimated  to  be  0.1  and  F  is  obtained  from 
the  Hertz-Knudsen  equation: 


F  = 


P-P 


(2mnkT) 


1/2 


(12) 


where  P  is  the  partial  pressure  of  zinc,  Pe  is  the  equilibrium  partial  pressure  of  zinc  at  the 
deposition  location,  m  is  the  mass  of  zinc  atom,  and  T  is  the  local  temperature  (1670  K  in  this 
case).  Because  there  was  no  deposition  prior  zinc  vapor  reached  the  growth  region  of  plate-caped 
micropillar  structure,  we  can  use  the  total  amount  yielded  by  the  precursor  as  the  initial  partial 
pressure  of  zinc  to  calculate  the  condensation  rate  of  zinc  for  the  formation  of  plate-caped 
micropillar  structure. 

The  unit  for  zinc  condensation  rate  Rcon  is  mol/(m2-s).  It  gives  the  deposition  amount  per 
second  per  square  meter.  For  an  infinitesimal  distance  dx,  the  deposition  time  is  dx/r,  where  r  is 
the  flow  rate  of  the  carrier  gas.  The  flow  rate  was  obtained  from  the  flowmeter  of  the  CVD 
system,  and  then  converted  to  the  flow  rate  in  the  growth  chamber  considering  the  pressure  and 
temperature  in  the  chamber  by  a  rough  estimation  using  PV  =  nRT.  The  deposition  area  is  the 
internal  surface  of  the  tubular  chamber,  n-D  ■ dx ,  where  D  is  the  internal  diameter  of  the  furnace 
tube.  The  deposition  volume  is  n{F>/2f dx.  Thus,  the  amount  of  deposition  at  this  infinitesimal 
area  per  unit  volume  (i.e.  the  vapor  concentration  reduction)  is 


7r-D-dx-  —  -R, 


Ac  = 


n  ■  (— )2  •  dx 
2 


Eq.  13  is  rearranged  and  simplified  as: 


Ac  = 


jKon 
D ■  r 


■  x 


(13) 


(14) 


where  x  is  the  distance  from  the  deposition  starting  point  to  the  point  of  interest.  The  deposition 
starting  point  can  be  measured  from  the  experiment. 


Using  P  =  cRT  to  convert  the  concentration  to  pressure,  we  can  combine  Eq.  14  and  a  = 
P/Pe  -  1  to  get: 


(7  = 


(co  ~ 


4PC„„ 

D  •  r 


■  x)RT 


—  1 


Arranging  Eq.  15,  we  have: 


cr  = 


c0RT~Pe 


A 

PeDr 


(15) 


(16) 
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Through  this  equation,  the  position  of  the  microplate  can  be  converted  to  the  local 
supersaturation. 


Figure  2.  Wedding  cake  growth  model.  A.  Microplates  with  various  numbers  of  terraces.  The 
number  of  terraces  decreased  as  the  sampling  point  moved  along  the  gas  flow  (downstream) 
direction,  i-ix  are  microplate  samples  collected  at  position  16.8479  cm,  16.8664  cm,  16.8669  cm, 
16.9033  cm,  16.9189  cm,  16.9256  cm,  16.9459  cm,  16.9741  cm,  and  17.0129  cm  from  the  center 
point,  respectively.  All  the  scale  bars  are  1  pm.  B.  Plots  of  the  axial  to  lateral  growth  rate  ratio  as 
a  function  of  supersaturation.  Black  dots  are  measured  data.  Black  line  is  the  fitted  curve  based 
on  the  measured  data,  where  a  linear  relationship  can  be  determined.  Blue  line  is  derived  from 
wedding  cake  growth  model.  Red  line  calculated  is  from  screw  dislocation  growth  model  for 
comparison. 

This  relation  shows  that  the  va/vi  ratio  increases  with  supersaturation,  suggesting  a  trend  that 
the  number  of  terrace  increases  with  supersaturation  in  Figure  2B.  The  significantly  closer  match 
of  the  wedding  cake  model  over  screw  dislocation  model  further  evidences  the  governing  role  of 
the  wedding  cake  mechanism  in  the  evolution  of  ZnO  microplate  structure. 

In  order  to  reveal  the  er-related  evolution  of  wedding  cake  tip  morphology  in  a  more  general 
NW  growth  scenario,  the  deposition  condition  was  adjusted  to  reach  desired  a  for  NW  formation. 
Under  this  condition,  ZnO  NWs  were  formed  with  uniform  size  and  geometry,  where  each  NW 
exhibited  a  terraced  tip  feature.  Just  like  the  microplates,  the  central  terrace  exhibited  a  flat 
surface  and  no  hint  of  the  existence  of  any  screw  dislocation  could  be  observed  by  both  SEM  and 
TEM  characterizations.  All  the  observations  evidenced  the  wedding  cake  growth  mechanism  in 
the  formation  of  these  NWs.  Furthermore,  the  tip  terrace  feature  exhibited  the  same  relation  with 
the  deposition  temperature  as  that  observed  from  microplates.  When  the  deposition  temperature 
decreases,  the  number  of  terrace  at  the  NW  tip  decreases  and  eventually  the  flat  tip  morphology 
was  obtained.  After  a  relative  large  temperature  range,  the  suspended  microplate  morphology 
appeared,  suggesting  the  deposition  condition  reached  the  critical  low  o  region.  Similar  change 
could  also  be  observed  by  increasing  the  oxygen  partial  pressure.  Qualitatively,  high  oxygen 
content  in  the  reaction  chamber  hinders  the  decomposition  of  ZnO  precursors  and  thus  lowers  a. 
Therefore,  as  the  oxygen  partial  pressure  increased,  the  NW  tip  exhibited  the  similar  change 
from  a  large  number  terraced  surface  to  a  flat  tip. 


2.  Observation  of  inverted  wedding  cake  growth  mechanism 

Another  wedding  cake-related  growth  phenomenon  was  observed,  which  was  similar  to  the 
step  edge  growth  in  VSS  processes.  A  series  of  single  crystal  hexagonal  branch-crowned  ZnO 
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NW  structure  was  obtained  by  dynamically  changing  the  deposition  temperature  during  the 
growth.  The  top  structure  was  found  to  be  a  basin-like  concave  crown  with  six  equivalent 
extended  branches.  An  inverted  wedding  cake  growth  model  involving  the  Ehrlich-Schwoebel 
(ES)  barrier  and  step  propagation  along  the  top  surface  edge  was  developed  to  explain  the 
growth  mechanism.  By  changing  the  oxygen  partial  pressure  and  deposition  temperature,  the  ES 
barrier  of  ZnO  was  quantified. 


Figure  3.  Morphologies  of  crown-caped  NW  structures.  A.  Large  scale  side  view  and  B.  top 

view  of  the  crown-caped  NW  arrays.  C.  A  tilted  top  view,  D.  a  side  view,  and  E.  a  tilted  bottom 
view  of  an  individual  six-branched  crown  cap  structure.  F.  A  top  view  of  the  crown  surface 
showing  the  hexagonal  center  pattern.  G.  Enlarged  image  on  the  crown  center,  showing  the 
hexagonal  concentric  steps.  H.  Enlarged  image  on  one  branch,  showing  a  large  flat  surface  area 
around  the  tip  of  the  branch  and  the  steps  became  more  regularly  spaced  near  the  center  region. 

Dynamically  controlled  VS  growth  of  ZnO  nanostructures  was  carried  out  under 
dynamically-tuning  deposition  conditions.  Both  the  evaporation  and  deposition  temperature  was 
increased  during  the  growth  providing  a  non-steady  state  growth  environment.  The  rapid  change 
of  vapor  concentration  by  the  end  of  growth  period  yielded  a  unique  morphology-switch  growth 
phenomenon,  where  all  the  NWs  were  capped  with  a  significantly  broadened  crown,  as  shown  by 
the  side-viewed  SEM  image  in  Figure  3A.  Different  from  our  previous  observations  of  the 
hexagon-capped  nanopillars,  all  the  top  crowns  exhibited  clear  six-fold  branches  with  uniform 
sizes  and  distributions  (Fig.3B).  Figure  3C  shows  a  typical  structure  of  the  hybrid  morphology. 
The  NW  had  a  diameter  of  -500  nm  and  a  length  of  -30  pm  or  even  longer.  The  crown  had  a 
core  diameter  of  -3  pm  and  six  branches,  which  expanded  up  to  -7  pm  wide.  Side  view  in  Fig. 
3D  shows  that  the  crown  had  a  basin-like  shape,  where  the  diameter  gradual  increased  from  -500 
nm  at  the  NW  junction  to  -7  pm  at  the  top.  When  observed  from  the  bottom  as  shown  in  Fig.  3E, 
the  junction  between  the  NW  and  crown  showed  a  very  smooth  transition  and  no  discontinuity 
could  be  identified. 
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The  top  view  in  Fig.  3F  reveals  more  details  about  the  crown’s  surface  features.  The  crown 
center  had  a  hexagonal  outline  as  marked  by  the  white  dashed  line  and  six  branches  grew  out 
from  the  six  edges  of  the  hexagon.  The  enlarged  view  of  the  plate  center  in  Fig.  3G  shows  the 
surface  was  featured  by  concentric  hexagonal  circular  steps.  At  the  innermost  circle,  there  was  a 
flat  surface  and  no  screw  dislocation  spiral  center  could  be  identified.  Different  from  our 
previous  observations,  the  step  height  increased  from  the  plate  center  to  plate  edge,  forming  a 
concave  surface  instead  of  the  classic  mount-shaped  feature.  Similar  step  structure  could  also  be 
observed  on  the  branch  surface  (Fig.  3H).  At  the  very  front  branch  tip  region,  a  relatively  wide 
flat  surface  area  was  commonly  observed,  which  could  be  more  than  1  pm  wide.  The  initial  steps 
exhibited  larger  step  width  than  the  steps  in  the  center  region.  The  front  of  each  step  was  not 
straight  and  does  not  represent  any  crystal  facets.  The  flat  area  also  covered  a  small  area  along 
the  edge  of  the  branches  as  marked  by  dashed  lines,  indicating  edge  nucleation  and  growth  was 
much  more  favorable  than  other  positions  on  the  plate.  These  morphological  features  suggested 
that  new  atomic  layers  would  nucleate  first  at  the  edge,  propagate  inward,  and  slow  down  when 
approaching  the  previous  step  due  to  the  limited  quantity  of  adatoms  on  the  surface. 

The  crystal  structure  of  the  crown  and  junction  region  was  characterized  by  transmission 
electron  microscopy  (TEM).  HRTEM  images  were  taken  from  the  joint  part  between  the  crown 
and  the  NW.  It  showed  that  the  lattice  was  continuous  at  the  connection  and  no  dislocation  core 
could  be  observed.  The  SAED  patternalso  confirmed  the  single  crystal  feature  of  the  entire 
hybrid  structure. 

In  this  model,  the  layer-by-layer  growth  was  initiated  at  the  plate  edge  instead  of  nucleating 
at  the  center  as  regular  wedding  cake  growth.  The  newly  formed  atomic  layers  then  propagated 
toward  the  plate  center  as  adatoms  adsorbed  and  diffused  toward  the  nucleation  sites,  which  thus 
formed  a  circular  step  along  the  plate  edge.  The  expansion  of  the  top  surface  is  a  result  of  the 
diffusion  down  of  the  adatoms  by  overcoming  the  ES  barrier.  Since  the  ES  barrier  prevents  the 
adatoms  from  overcoming  the  crystal  ledge,  it  helps  increase  the  concentration  of  surface 
adatoms  and  thus  facilitates  the  nucleation  of  new  atomic  layers.  According  to  our  previous 
results,  as  the  supersaturation  (a)  decreases,  the  chance  for  adatoms  to  overcome  the  ES  barrier 
rises,  while  the  concentration  of  surface  adatoms  drops.  At  a  certain  small  a,  the  top  surface 
would  become  no  longer  a  favorable  place  for  nucleation  of  new  atomic  layers  due  to  the 
reduced  adatom  concentration.  Instead,  the  adatom  concentration  rises  at  the  crystal  edge  which 
becomes  a  new  nucleation  site.  Therefore,  new  atomic  layers  start  to  form  along  the  edge  of  the 
top  surface  while  it  becomes  wider  as  the  adatoms  continuously  attach  to  the  edge.  In  typical 
mount  morphology,  the  top  layer  has  the  largest  flat  area  due  to  the  slower  nucleation  rate 
compared  to  the  ledge  expansion  rate  (nucleation-limited  region).  Similarly,  in  the  inverted 
wedding  cake  growth,  the  top  most  layer,  which  is  the  surface  along  the  edges  of  the  branches, 
also  exhibited  the  largest  flat  surface  area  (Fig.  2H).  As  the  steps  approaching  the  center  area,  the 
space  between  them  became  much  smaller.  This  is  consistent  with  the  fact  that  the  step  growth 
slows  down  when  it  approaches  the  growth  front  of  its  previous  one  due  to  the  decrease  of  the 
effective  adatom  diffusion  area,  suggesting  a  mass-transport  limited  growth  region.  The  large 
step  height  also  suggested  that  the  step  bunching  effect  was  a  common  growth  feature  in  ZnO 
layer-by-layer  growth  due  to  the  high  possibility  for  adatoms  to  overcome  the  ES  barrier  during 
steps  propagation. 

In  order  to  further  understand  how  the  inverted  wedding  cake  growth  is  related  to 
conventional  wedding  cake  growth  of  NWs,  a  series  of  syntheses  with  different  oxygen  partial 
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pressures  were  conducted.  The  oxygen  partial  pressures  and  the  temperature  zones  would  modify 
both  the  condensation  and  oxidation  rates  of  zinc,  and  thus  influence  the  growth  kinetics.  In  the 
highest  temperature  region,  the  nanostructures  appeared  in  the  form  of  NWs  with  a  clear  mount¬ 
shaped  top  surface.  In  the  medium  temperature  zone,  a  hexagonal  plate  appeared  at  the  top  of 
each  NW  and  a  concave  center  started  to  develop.  As  temperature  further  reduced,  the  top  plate 
clearly  showed  the  basin  shape  and  six  equivalent  branches  evolved  and  elongated  from  the  plate 
edge.  The  size  of  the  plates  also  increased  when  lowering  the  deposition  temperature.  When 
applying  different  oxygen  partial  pressures,  the  same  trend  was  observed.  NW  morphology  was 
obtained  at  the  lowest  oxygen  partial  pressure.  Increasing  the  oxygen  partial  pressure  led  the  top 
morphology  evolution  from  mount-shape  hexagonal  plates  to  basin-shaped  branched  structures. 

The  unique  morphology  switching  phenomenon  could  be  understood  by  the  dynamic 
deposition  condition  change  during  the  growth  processes.  At  the  beginning,  the  deposition 
temperature  was  relatively  low,  which  keeps  the  equilibrium  partial  pressure  (Pe)  at  a  relatively 
low  level.  Since  the  precursor  vapor  concentration  ( P )  remained  at  a  relatively  high  level,  a  at 
any  deposition  temperature  was  large.  This  condition  favored  NWs  growth  through  the  layer-by- 
layer  wedding  cake  model  along  the  [0001]  direction.  At  this  stage,  the  adatoms  were  confined  at 
the  top  (0001)  surface  and  the  diffusing  down  was  significantly  restricted  by  the  ES  barrier.  As 
the  deposition  temperature  rises,  Pt  increased  rapidly  leading  to  a  quick  drop  of  o.  The  reduction 
of  a  was  the  most  prominent  at  the  highest  deposition  region,  leaving  insufficient  time  for  crown 
formation.  Thus,  NWs  with  a  mount- shaped  tip  was  obtained  as  the  final  product.  Meanwhile, 
relatively  small  a  was  still  remained  at  lower  deposition  temperature  zones.  Such  a  condition 
favored  the  diffusion  down  of  the  adatoms  from  the  top  surface  by  overcoming  the  ES  barrier 
and  nucleated  an  additional  layer  circling  around  the  tip  of  the  NW  broadening  the  top  layer.  As 
the  diffusion  down  proceeded,  the  top  surface  became  undersaturated  for  nucleating  new  atomic 
layers  and  the  growth  was  completely  initiated  from  the  edge.  Thus,  the  growth  mode  switched 
to  the  inverted  wedding  cake  growth  and  the  hexagonal  crowns  were  formed.  The  growth  could 
be  terminated  at  this  stage  when  Pe  further  increased  to  diminish  a  at  the  medium  deposition 
temperature  zone.  However,  at  the  even  lower  deposition  temperature  region,  a  could  be 
favorable  for  lateral  growth  and  larger  crowns  were  obtained  in  that  region.  The  appearance  of 
six  branches  might  be  a  result  of  oxygen  adsorption  that  selectively  restricted  the  growth  on  the 
flat  edge  area.  High  oxygen  partial  pressure  hinders  the  decomposition  of  ZnO  precursor  and 
thus  lowers  a.  Hence,  the  high  oxygen  partial  pressure  favors  the  growth  of  the  crown  featured 
morphologies  and  also  promotes  the  formation  of  branched  structure. 
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Figure  4.  Relationship  between  the  growth  rate  of  microplates  and  temperature  under  different 
oxygen  partial  pressure,  1  Pa  (black  dots),  10  Pa  (red  dots),  30  Pa  (blue  dots)  and  50  Pa  (green 
dots).  The  correspondingly  colored  lines  are  fitted  curves  which  give  similar  slopes  revealing  the 
activation  energy  for  overcoming  the  ES  barrier.  Each  dot  is  the  average  value  of  ten  sets  of 
experimental  data  from  the  same  temperature  region  and  error  bars  show  standard  error. 

In  general,  the  plate  expansion  is  a  kinetic  process,  which  is  determined  by  the  diffusion 
down  rate  (a).  The  relationship  between  a  and  the  ES  barrier  ( EEs )  follows  the  equation: 

K  —  E 

a  =  ex p( — — - — )  (17) 

kT 

where  Esj  is  the  surface  diffusion  barrier,  k  is  the  Boltzmann  constant,  and  T  is  temperature. 
Therefore,  the  deposition  temperature-related  plate  morphology  can  be  used  to  identify  the  ES 
barrier  for  ZnO  nanostructure  growth.  From  Eq.  17,  the  growth  rate  of  the  suspended  plate  due  to 
adatom  diffusing  down  to  the  edge  (vj)  can  be  derived  as  a  function  of  T. 

lnvrf=  In  a  _(££s-^+M/  +9087j)l  (18) 

8^27mk  k  T 

where  a  is  the  distance  between  adsorption  sites,  AU  is  the  kinetic  barrier  for  the  incorporation  of 
the  building  units  into  the  half  crystal,  8  is  average  spacing  of  the  kink  sites,  and  m  is  the  mass  of 
zinc.  Vd  can  be  determined  by  comparing  the  suspended  length  and  the  corresponding  height  of 
plates,  which  represents  the  ratio  between  vertical  growth  rate  and  diffusion  down  rate.  Through 
this  method,  Vd  was  identified  from  each  temperature  zone  from  samples  grown  under  different 
oxygen  partial  pressures.  The  Arrhenius  plots  of  In  Vd  vs.  1/T  under  the  four  selected  oxygen 
partial  pressures  showed  that  the  linear  fits  of  the  data  exhibited  very  close  slopes,  from  which 
the  ES  barrier  EEs  was  found  to  be  -0.88  ±  0.33  eV,  in  Figure  4.  Compared  to  the  reported  values 
in  other  materials  (0.52  eV  for  Au,  0.40  eV  for  Cu  and  0.12  eV  for  Al),  the  value  we  identified 
was  much  higher.  This  could  be  the  reason  that  the  diffusion  down  phenomenon  was  not 
commonly  observed  in  oxide  nanocrystals. 


3.  Kinetics-Driven  Crystal-Facets  Evolution  at  the  Tip  of  Nanowires  -  A  New 
Implementation  of  the  Ostwald-Lussac  Law 

Through  controlling  the  growth  kinetics,  the  facets  evolution  involving  the  appearance  of  the 
high  index  facets  at  the  nanowire  tips  was  observed.  The  facets  evolution  was  induced  by  the 
change  of  the  relative  growth  rates  of  different  facets  and  it  followed  the  Ostwald-Lussac  Law. 
The  results  provide  a  platform  to  understand  the  nanowire  growth  mechanism  from  the  vapor 
phase. 
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Figure  5.  SEM  images  and  corresponding  schematic  showing  facets  evolution  at  NW  tips. 

(A)  -  (C)  NW  with  the  tip  enclosed  by  six  a  facets  and  one  facet  (0001).  From  (A)  to  (C),  facet 
(0001)  gradually  decreased  when  facets  a  increased  simultaneously.  (D)  -  (F)  NW  with  the  tip 
covered  by  six  a  facets  and  six  (1  facets.  From  (D)  to  (F),  facets  P  appeared,  expanded  and 
eventually  covered  the  entire  tip  area  at  the  expenses  of  facets  a.  (G)  -  (I)  NW  with  six  p  facets 
and  one  (0001)  facet  co-existed  tip.  From  (G)  to  (I),  facets  p  gradually  vanished  and  facet  (0001) 
appeared  again  with  a  30°  rotation  relative  to  the  six  sidewall  of  the  NW.  In  the  schematic, 
yellow  represents  facet  (0001),  orange  is  facets  a,  and  blue  is  facets  p.  All  scale  bars  are  100  nm. 


ZnO  NWs  were  grown  through  a  dynamically  controlled  vapor-solid  deposition.  During  the 
NW  growth,  both  the  evaporation  temperature  and  the  deposition  temperature  were  steadily 
increased  in  order  to  achieve  an  accelerated  dynamic  growth  environment  to  make  possible  the 
appearance  of  the  thermodynamically  unstable  morphology  growth.  Under  such  a  non-steady 
state  deposition  condition,  randomly-oriented  submicron-sized  ZnO  NWs  were  grown  on  the 
poly  crystalline  alumina  substrate  within  a  broad  temperature  range  from  1300  K  to  1673  K.  The 
NWs  were  typically  20  pm  long  with  a  diameter  of  -300  nm.  Most  NWs  exhibited  a  clear 
hexagonal  cross  section  due  to  the  P63W.C  lattice  symmetry  of  the  wurtzite  crystal  structure. 
Distinguishing  from  common  ZnO  NWs  that  typically  have  a  flat  hexagonal  (0001)  tip  facet,  the 
tips  of  ZnO  NWs  received  from  the  dynamic  deposition  condition  were  mostly  enclosed  by  at 
least  two  groups  of  crystal  facets.  These  tip  facets  also  distributed  with  a  hexagonal  symmetry. 
The  NW  tip  facets  evolution  were  only  found  at  a  high  deposition  temperature  region  (1665K) 
and  directly  related  to  the  deposition  location.  As  shown  by  the  relationship  between  deposition 
temperature  and  location  in  the  deposition  chamber  in  Figure  S2,  the  temperature  was  almost 
constant  within  that  region.  At  other  deposition  region,  the  NWs  tip  was  the  regular  flat  (0001) 
facet.  Within  this  narrow  deposition  region  window  (-1  cm),  different  crystal  facets  evolved  and 
vanished  as  the  deposition  location  changed  slightly. 

Representative  tip  morphology  illustrating  the  tip  facets  evolution  as  a  function  of  location 
and  corresponding  schematic  drawing  are  presented  in  Figure  5.  At  the  position  closest  to  the 
chamber  center  (17.7  cm  from  the  chamber  center,  where  the  precursor  source  was  located),  the 
NW’s  tip  was  only  partially  covered  by  the  hexagonal  (0001)  facet  at  the  center,  with  six  angled 
facets  (named  as  a  facet)  evenly  distributed  along  the  six  edges  between  the  vertical  sidewall 
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facet  and  the  top  (0001)  facet,  as  shown  in  Figure  5  A.  As  the  deposition  position  moved  slightly 
downstream  by  ~3  mm,  the  six  a  facets  gradually  expanded  and  eventually  occupied  the  entire 
tip  area.  Meanwhile,  the  (0001)  facet  shrank  into  a  sharp  point  at  the  center  of  the  tip  (Figure  5  A 
to  C).  Subsequently,  another  group  of  facets  (named  as  P  facet)  emerged  at  the  six  radial  comers 
of  the  tips  (Figure  5D).  The  p  facets  expanded  and  eventually  covered  the  entire  tip  area  at  the 
expenses  of  the  a  facets  as  the  deposition  position  moved  downstream  by  another  ~3  mm  (Figure 
5D  to  F).  At  the  position  slightly  farther  away  from  the  chamber  center,  the  hexagonal  (0001) 
facet  appeared  again  from  the  center  of  the  tip.  Different  from  the  original  (0001)  facet  obtained 
from  the  high-temperature  NWs,  however,  the  re-emerged  (0001)  facet  exhibited  a  30°  rotation 
relative  to  the  hexagonally-distributed  NW  sidewalls  (Figure  5G).  The  (0001)  facet  expended 
and  eventually  took  over  the  entire  NW  tip  again  at  the  deposition  position  of  18.7  cm  away  from 
the  chamber  center,  where  the  p  facets  vanished  completely  (Figure  5H,  I).  Across  this  narrow 
deposition  region,  the  deposition  temperature  difference  was  negligible  (<5°C).  The  vapor 
pressure  of  the  precursor  changed  significantly  due  to  the  rapid  deposition  across  this  region. 
These  suggested  that  the  precursor  supersaturation  was  the  dominating  factor  at  this  high 
deposition  temperature  area.  Thus,  the  facet  evolution  is  primarily  correlated  to  the  location- 
induced  supersaturation  difference. 


Figure  6.  SEM  images  and  TEM  projections  for  facets  indexing.  (A)  NW  with  facets  {112  2} 
and  facets  {1013}  co-existed  tip.  TEM  projection  and  selected  area  electron  diffraction  (SAED) 
patterns  (inset)  along  (B)  the  direction  [21 10]  and  (C)  the  direction  [1010],  (D)  NW  with  facets 
{1013}  and  facet  (0001)  co-existed  tip  and  (E)  the  corresponding  TEM  projection.  Inset  is  the 
SAED  pattern  along  [2110].  (F)  NW  with  tip  totally  solely  covered  by  facets  {1122}  and  (G) 
corresponding  TEM  projection.  Inset  is  the  SAED  pattern  along  [1010],  (H)  NW  with  facets  {11 
2  2}  and  facet  (0001)  co-existed  tip  and  (I-K)  corresponding  TEM  projection.  Insets  are  the 
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SAED  patterns  along  [1010]  showing  the  facets  {112  2}  and  facet  (0001)  evolution  under 
different  supersaturation. 

In  order  to  understand  the  tip  facets  evolution  phenomenon,  the  involving  facets  were 
indexed  first.  Transmission  electron  microscopy  (TEM)  was  implemented  to  characterize 
selected  NW  tips  from  various  projection  directions.  A  NW  with  both  the  a  and  [3  facets  was 
imaged  as  shown  in  Figure  6A.  TEM  image  was  first  recorded  along  the  [2110]  projection 
direction,  where  only  the  a  facets  can  be  seen  (Figure  6B).  The  projection  plane  was  also 

outlined  by  the  dashed  blue  line  in  Figure  6 A.  Thus,  the  angle  between  the  {1010}  and  a  facets 
was  directly  measured  to  be  121°  from  the  projected  image,  which  matches  well  to  the  selective 
area  electron  diffraction  (SAED)  pattern  (inset  of  Figure  6B).  It  is  noted  that  in  the  SAED  pattern, 
presence  of  the  forbidden  Bragg  spots  of  the  (0001)  and  (0003)  planes  was  attributed  to  the 

secondary  Bragg  electron  diffraction  according  to  the  satisfied  condition  ( e.g .:  (01 11)  +  (01 10)  = 
(0001)).  Therefore,  the  a  facet  was  indexed  to  be  facets  {1013},  a  high  index  plane.  By  tilting 
the  same  NW  by  30°  along  the  NW’s  axial  axis  (the  [0001]  direction),  another  projection 
consisting  of  both  the  a  and  p  facets  were  imaged,  as  shown  in  Figure  6C.  The  SAED  pattern 

(the  inset  of  Figure  6C)  confirmed  that  projection  was  now  along  the  [1010]  direction.  The 
geometry  of  the  projected  image  was  outlined  by  the  dashed  red  line  in  Figure  6 A.  At  this 
condition,  the  P  facets  were  parallel  to  the  electron  beam,  enabling  the  measurement  of  the  angle 
between  the  p  facets  and  the  horizontal  plane  (known  to  be  the  (0001)  plane).  The  angle  was 
measured  to  be  122°.  Therefore,  from  the  corresponding  SAED  pattern  (inset  of  Figure  6C),  the 

P  facet  was  determined  to  be  {112  2},  another  high  index  plane. 

NWs  at  different  evolution  stages  were  collected  to  further  confirm  the  combination  of  tip 
facets.  One  NW  with  the  {1013}  facets  co-existed  with  the  (0001)  facet  at  the  tip  was  imaged 

and  characterized  (Figure  6D).  A  TEM  image  was  taken  from  a  projection  plane  of  {2110} 
(Figure  6E),  as  outlined  by  the  dashed  blue  line  in  Figure  6D.  The  orientation  was  validated  by 
the  corresponding  SAED  pattern  (inset  of  Figure  6E),  which  was  the  same  as  the  SAED  pattern 

of  the  projection  in  Figure  6B.  The  angle  between  the  {1013}  facet  and  the  {1010}  facets  was 
measured  to  be  121°,  identical  to  that  value  measured  in  Figure  6B,  confirming  the  same  facets 

configuration  in  this  NW  tip.  Another  NW  with  the  tip  solely  covered  by  six  {112  2}  facets 

(Figure  6F)  was  imaged  along  the  [1010]  direction  (Figure  6G).  The  angle  between  the  [112  2] 
facet  and  the  horizontal  plane  (dashed  line,  (0001))  was  measured  to  be  122°,  indicating  the  same 

facet  group  as  those  shown  in  Figure  6C.  Another  NW  with  the  combination  of  [112  2]  and 
(0001)  facets  at  the  tip  is  shown  in  Figure  6H.  TEM  image  projected  along  the  [1100]  direction 
is  shown  in  Figure  61.  The  angle  between  the  flat  tip  and  the  {112  2}  facet  was  measured  to  be 
122°,  identical  to  the  measurement  when  the  (0001)  facet  was  not  adjacent  (Figure  6C).  All  these 
measurements  confirmed  that  there  were  the  same  groups  of  crystal  facets  during  the  tip 
morphology  evolution.  The  TEM  image  also  showed  clear  thickness  fringes  around  the  NW 
edges.  The  different  fringe  patterns  from  outside  to  the  center  directly  correlated  to  the  facets 
organization.  Other  NWs  obtained  at  slightly  farther  away  from  the  chamber  center  that  had  the 
same  facets  combination  but  increasing  (0001)  areas  are  also  shown  in  Figure  6J  and  6K.  All  the 
SAED  patterns  and  facets  angles  were  identical,  confirming  no  more  new  facets  were  evolved 
when  the  NW  tip  was  reclaimed  by  the  (0001)  facet  again. 
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Figure  7.  High  resolution  TEM  images  on  the  lattice  structure  at  the  NW  tip.  (A)  NW  with 
the  tip  predominately  covered  by  facets  {1013}  and  (B)  the  corresponding  TEM  projection.  Inset 
is  the  SAED  pattern  along  [303  1].  (C)  HRTEM  at  the  tip  region,  marked  in  (B),  showing  facets 

{1013}  sharp  connected  by  a  small  portion  of  facet  (0001).  (D)  Facet  {1013},  marked  in  (B), 
showing  the  rough  surface.  Islands  with  a  few  atomic  layer  thickness  covered  almost  the  entire 

surface.  (E)  Corner  between  facet  {1013}  and  {1010},  marked  in  (B),  showing  multi-atomic- 
layer  steps,  induced  by  the  surface  energy  minimization. 

High  resolution  TEM  (HRTEM)  images  were  obtained  from  the  NW  tip  (shown  in  Figure 
7 A)  to  reveal  the  lattice  structure  around  the  tip  edges.  The  selected  NW  had  a  tip  predominately 

covered  by  the  {1013}  facets  (Figure  7B).  In  general,  the  NW  exhibited  a  perfect  crystal  lattice 
without  any  observable  dislocations.  No  hint  of  screw  dislocation  could  be  identified.  From  the 
very  center  area,  as  shown  in  Figure  7C,  a  small  portion  (~3  nm)  of  the  (0001)  facet  could  still  be 

observed.  It  quickly  turned  into  the  {1013}  facets  with  sharp  comers.  The  lattice  spacing  of  the 
{1013}  plane  was  measured  to  be  0.1425  nm,  consistent  with  0.1477  nm  calculated  from  the 

JCPDS  card  (No.  36-1451).  The  {1013}  facets  were  found  to  be  rough,  as  shown  in  Figure  7D. 
Islands  with  a  few  atomic  layer  thickness  and  multiple  atomic  steps  covered  almost  the  entire 
surface.  This  could  be  a  result  of  the  surface  energy  minimization  of  the  crystal  system.  For  the 

same  reason,  the  {1013}  facet  restructured  into  a  series  of  multi-atomic-layer  steps  at  the 

junction  with  the  {1010}  facet.  The  rough  surface  and  corners  might  be  responsible  for  the 
crystal  facet  evolution  during  the  NW  growth  at  different  temperature  and  supersaturation,  i.e. 
different  driving  force  of  atom  diffusion. 

Above  crystal  lattice  analysis  suggested  that  the  NWs  were  likely  grown  by  a  layer-by-layer 
mechanism,  as  we  discovered  previously.  In  this  mechanism,  the  Ehrlich-Schwoebel  (ES)  barrier 
is  a  critical  parameter  that  prevents  the  adatoms  from  diffusing  over  the  atomic  steps  and  thus 
builds  up  necessary  adatom  concentration  to  nucleate  new  atomic  layers.  In  an  ideal  ZnO  NW 
growth  process,  atomic  layers  of  the  (0001)  plane  nucleate  at  the  NW  tip  and  grow  laterally 
covering  the  entire  tip  surface.  Nevertheless,  as  the  (0001)  plane  approaches  the  edge,  e.g.  the 
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comer  with  the  (1010)  plane,  new  high-index  crystal  facets  will  form,  such  as  the  {1013}  and 

{112  2}  families  (Figure  8 A).  These  facets  have  high  surface  energy  and  typically  vanish  rapidly 
as  their  sizes  become  appreciable  (e.g.  a  few  nm  wide).  Thus,  these  high-index  planes  can  be 
considered  as  metastable  intermedia  stages  as  the  complete  evolvement  of  a  new  (0001)  plane. 
However,  our  observation  of  the  co-existence  of  these  high-index  facets  within  a  narrow 
deposition  window  suggested  that  under  certain  circumstances,  these  intermedia  stages  might 
stably  present.  For  example,  the  ES  barrier  for  adatoms  to  diffuse  from  the  (0001)  plane  to  the 

(1013)  plane  might  be  significantly  higher  than  that  for  adatoms  to  diffuse  backward.  Therefore, 
the  overall  energy  barrier  for  the  formation  of  a  complete  new  (0001)  plane  would  be  much 

higher  than  that  for  leaving  the  high-energy  (1013)  plane  exposed.  Similar  scenario  would  also 

apply  to  the  facets  transition  between  the  pairs  of  (000 1 )/ { 1 1  2  2}  and  {1013}/{1122}  as  we 
observed  experimentally. 


Figure  8.  Analogy  of  facet  evolution  to  the  Ostwald-Lussac  Law.  (A)  Schematic  of  the  facet 
formation  model.  (B)  Schematic  energy  landscape  of  a  series  of  facets  illustrating  the  facets 
evolution  pathway.  At  certain  supersaturation,  the  high-index  facets  tend  to  be  stable  due  to  the 
lower  energy  barriers.  (C)  Dependence  of  the  growth  rate  ratio  (Left  y-axis)  and  the  energy 
barrier  differences  (Right  y-axis)  of  (0001)/{  1013},  { 101 3}/{1122}  and  (0001)/{  11 2 2}  on  the 
supersaturation. 

This  free  energy  barrier-induced  appearance  of  metastable  intermediate  stages  shares  much 
similarity  with  the  Ostwald-Lussac  (OL)  Law  that  has  been  primarily  discovered  in 
biomineralization.  The  OL  law  is  an  empirical  law  of  phase  transformation  stating  that  an 
unstable  system  does  not  necessarily  transform  directly  into  the  most  stable  phase,  but  likely  into 
one  with  the  minimum  energy  loss.  Here,  we  propose  that  the  crystal  facet  formation  could  also 
be  represented  by  a  similar  energy  landscape  of  the  OL  law,  although  the  atomic  plane  growth  is 
not  really  a  phase  transformation.  In  analogy  to  the  OL  law,  a  free  energy  landscape  was 
constructed  to  show  the  ZnO  NW  tip  facets  evolution  relationships  based  on  our  experimental 
observations. 
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As  shown  in  Figure  8B,  stage  1  to  stage  4  represent  the  process  of  a  complete  (0001)  atomic 
layer  formation.  Stage  4  has  the  lowest  free  energy,  because  of  the  lowest  surface  energy  of  the 
(0001)  plane.  In  between  stages  1  and  4,  another  two  intermediate  stages  exist,  e.g.  formation  of 

the  {1013}  facets  (stage  2)  and  formation  of  the  {112  2}  facets  (stage  3).  These  stages  have  high 
energy  than  the  final  stage  4  because  of  the  higher  surface  energy  of  the  high-index  facets. 
Kinetically,  the  corresponding  energy  barriers  between  each  stage  were  marked  as  AGab,  where  a, 
b  represent  the  initial  and  final  stage,  respectively  during  the  crystal  facet  growth.  As  discussed 
previously,  these  energy  barrier  heights  might  change  toward  different  directions  as  the  ES 
barrier  varies  as  a  function  of  supersaturation.  Under  certain  circumstances,  the  energy  barriers 
involving  intermediate  stages  might  be  significantly  lower  than  AG14,  and  thereby  the  high-index 
facets  appear  in  the  final  NW  morphology.  The  relative  value  between  different  AG  dictates  the 
appearing  crystal  facet(s).  To  more  specifically  elaborate  this  empirical  law  of  facet  evolution, 
all  the  experimentally  observed  NW  tip  morphology  was  correlated  to  AG  relationships  and 
summarized  in  Table  1.  This  table  could  be  further  enriched  if  more  new  NW  tip  facets  are 
discovered. 

Similar  to  the  OL  law  for  phase  transformation,  our  proposed  law  for  facet  evolution  at  the 
NW  growth  front  is  purely  qualitative.  Quantification  of  all  the  relative  energy  barriers  requires  a 
comprehensive  knowledge  of  various  associated  energies  (e.g.  the  ES  barrier, 
adsorption/desorption  energy,  surface  energy,  etc.)  and  other  parameters  (e.g.  supersaturation, 
vibration  frequency,  etc.),  which  are  beyond  the  capacity  of  this  work.  Fortunately,  the  facet 
areas  demonstrated  a  good  relation  with  the  supersaturation.  According  to  kinetic  Wulff  s 
theorem,  the  relative  growth  rates  of  the  facets  can  be  deduced  from  the  relative  sizes  of  them. 
As  discussed  earlier,  the  supersaturation  (a)  was  considered  as  the  dominating  variable  in 
correlation  to  the  deposition  location  within  this  small  window.  Therefore,  based  on  the  evolved 

sizes  of  the  (0001),  {1013}  and  {112  2}  facets  measured  at  different  deposition  locations,  the 
relative  growth  rate  was  calculated  and  plotted  as  a  function  of  <7  (given  by  the  left  y-axis  in 
Figure  8C).  As  the  supersaturation  decreased,  the  relative  growth  rate  ratio  of  (0001)/{1013} 
decreased  from  1 . 02  to  0 . 2 1 ,  and  then  the  relative  growth  rate  ratio  of  { 101 3}/{ll  22}  decreased 

from  2.38  to  0.24,  and  finally  the  relative  growth  rate  ratio  of  (000 1)/{1122}  increased  from 
0.45  to  31.92.  The  relative  growth  rate  ratio  changes  indicated  that  the  growth  energy  barrier 
difference  between  the  corresponding  two  groups  of  facets  changed  at  different  supersaturation. 

According  to  the  crystal  growth  theory,  the  facet  growth  is  linearly  related  to  the 
supersaturation,  and  exponentially  related  to  the  temperature  and  the  growth  energy  barrier  of 
corresponding  facets,  which  can  be  expressed  by  the  equation: 

E 

v^Acrex p(--=)  (1) 

kl 

where  vi  is  the  growth  rate  of  facet  1,  A  is  a  parameter  related  to  the  growth  condition,  Ej  is 
the  growth  energy  barrier  of  the  facet  1,  k  is  Boltzmann  constant  and  T  is  the  deposition 
temperature.  Using  this  equation,  the  energy  barrier  differences  were  quantified  as  a  function  of 
a  (given  by  the  right  y-axis  of  Figure  8C).  For  all  of  these  three  pairs  of  facets,  the  energy  barrier 
difference  might  be  fitted  linearly  to  the  logarithm  of  the  supersaturation,  indicating  the  energy 
barrier  difference  was  logarithmically  influenced  by  the  supersaturation.  At  higher 
supersaturation  level,  corresponding  to  the  location  closest  to  the  chamber  center,  the  energy 
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barrier  for  the  formation  of  the  (0001)  facet  was  almost  the  same  as  that  of  the  {1013}  facets 
(only  1  meV  difference).  With  the  supersaturation  decreased,  the  energy  barrier  of  forming  the 

(0001)  facet  became  much  larger  than  that  of  the  {1013}  facets,  confirming  the  rapid  evolution 
of  the  {1013}  facets  in  place  of  the  (0001)  facet.  Similarly,  when  the  supersaturation  continued 
to  decrease,  the  energy  barrier  for  facets  {1013}  changed  from  being  50  meV  smaller  than  that 
of  the  {112  2}  facets  to  being  90  meV  larger  than  that  of  the  {1122}  facets,  providing  a  smaller 

energy  barrier  for  the  facet  evolution  from  {1013}  to  { 1 1  2  2 }.  As  the  supersaturation  decreased 
further,  the  energy  barrier  for  (0001)  formation  rapidly  decreased  from  50  meV  to  -210  meV 

relative  to  the  energy  barrier  for  {112  2}  formation,  confirming  a  relatively  larger  drive  force  for 
the  (0001)  to  quickly  cover  the  entire  NW  tip,  forming  a  flat  tip  surface. 
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